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Abstract 
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1 Introduction 



The problems of random surface growth have been an important subject of non-equilibrium 
statistical physics [TP] . Interestingly, surfaces show universal statistical behaviors, depend- 
ing on which mechanism plays a major role in growth. Among them is the Kardar-Parisi- 
Zhang (KPZ) universality class 3J. Though in many cases it is not enough to describe 
the growth of real materials, the KPZ universality class gives a satisfactory understanding 
of many growth models such as the Eden model, and hence plays a prominent role in 
theoretical study of growing surfaces. 

It is difficult to study the KPZ universality class analytically for general dimension. In 
one spatial dimension, however, some exact results had been obtained. For instance, the 
exponents are already obtained in [3| based on dynamical renormalization group techniques. 
These have been further confirmed by exact solutions for the models in the KPZ universality 
class . We will also restrict our attention to one spatial dimensional case in this article. 

Recently, more refined information about the fluctuation properties of the 1+1 dimen- 
sional KPZ universality class have been obtained 131111313110111111121131111113 • For some 
important quantities such as the height fluctuation, not only the exponents but also the 
scaling functions are obtained. The whole new developments originate from the work jTHj 
on the statistics of the longest increasing subsequence in a random permutation, which has 
turned out to be closely related to the random surface growth models [T7] . 

The polynuclear growth (PNG) models are known to belong to the KPZ universality 
class The PNG models are simple models which describe layer by layer random 

surface growth. In its standard version, a nucleation of a layer with height one occurs with 
rate two per unit length. After a nucleation, the layer grows laterally in both directions 
from the nucleation point with unit speed. In addition, there occur nucleations with rate 
two on already existing layers. A discrete version of the PNG model was introduced in JTSj 
and studied further in J3- The discrete PNG model has a close relationship with the 
asymmetric simple exclusion process (ASEP) [H|.ll3j . 

In P 1| 7 1| H] , the height fluctuation of the PNG model at a given time and a given position 
was studied. It turned out that the fluctuation strongly depends on the geometry of the 
models. For instance, for the model in the infinite space, the fluctuation at the origin is 
described by the GUE Tracy- Widom distribution in random matrix theory [2*312 1| . On the 
other hand, for the model in a half space with an external source at the origin, it is described 
by the GOE/GSE Tracy- Widom distributions 22 j or by the Gaussian. Models with other 
geometries are also considered in p"JIH|- These are based on the results in (3123112111213 > i n 
which the longest increasing subsequence problem with some symmetries are considered. 

For the models in the infinite space, multi-point equal time height fluctuations are 
studied by introducing the multi-layer version of the PNG model in [H] for the continuous 
time case and in [JJ"] for the discrete time case. It is shown that the scaling limit is described 
by the Airy process, which is closely related to the dynamics of the rightmost particle of 
the Dyson's Brownian motion model [23- 

In this article we study the multi-point equal time height fluctuations of the PNG model 
in a half space with external source at the origin. As mentioned above, the fluctuation 
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at a single point has already been considered in [23IZ1- At the origin, it is given by the 
GOE/GSE Tracy- Widom distribution whereas at other points, it is the GUE Tracy- Widom 
distribution. But the multi-point fluctuation has not been studied. Following the ideas 
of |27|I14|. we introduce the multi- layer PNG model. In the scaling limit we will show that 
the crossover between the fluctuation at the origin and that in the bulk corresponds to the 
orthogonal/symplectic to unitary transition at soft edge in random matrix theory [25] , 

The outline of the paper is as follows. In the next section, the model we consider is 
defined. In section |3J a determinantal process is considered. Section 0] deals with the 
scaling limit of the models. In section the results for the model without external source 
at the origin are summarized. The standard PNG model in a half space can be considered 
by taking the limit for the discrete model. This is done in section |3 Some discussions, 
conjectures and Monte-Carlo results are given in section [71 followed by the conclusion in 
the last section. 

2 Model 

Let r G N = {0, 1, 2, • • • } and t G N denote the discrete space and time coordinates 
respectively and h(r,t) the height of the surface at position r and at time t. Our model 
lives in a half space r G N, but we can extend h(r, t) to the whole space by setting 
h(r, t) = h(—r, t) for r G Z, i.e. the height is symmetric under the reflection with respect to 
the origin. We can further extend h(r, t) to all r G R by setting h(r, t) = h([r],t). Then the 
height looks symmetric under the reflection with respect to r = 1/2, i.e. h(r, t) = h(l — r, t) 
for r G R \ Z. The model considered in this article is described as follows. Initially at time 
t = we have a flat substrate. At an odd time t = 1, 3, • • • , nucleations occur at even sites 
r = 0, ±2, • • • , ±£ — 1. At the origin, a nucleation of a height k (G N) happens with the 
probability (1 — l^/q)iri^) k where < q < 1 and < 7 < l/^/q. For the height to be 
symmetric, nucleations at other points should be symmetric with respect to the origin; if 
the nucleation of a height k occurs at r, so does at — r as well. Each independent nucleation 
with height k occurs with probability (1 — q)q k . At an even time t = 2, 4, cdots, nucleations 
occur at odd sites r = ±1, ±3, • • • ± t — 1. The nucleations are again synchronous with 
respect to the origin and a height of each independent nucleation is a geometric random 
variable with parameter q. In the meantime the steps grow laterally in both directions with 
unit speed. Sometimes a downstep and an upstep collide with each other, in which case 
the two steps merge to one with the higher height. See Fig. 1. As is clear, the parameter 
q is related to the frequency of nucleations in the bulk whereas the parameter 7 represents 
the strength of the external source at the origin. The bigger 7 is, the stronger the source 
is. In particular 7 = corresponds to the model without the external source at the origin. 
Mathematically, our discrete PNG model can be defined by 

h{r, t + l) = max(/i(r - 1, t), h(r, t), h(r + 1, t)) + u(r, t + 1), (2.1) 

with the initial condition h(r, 0) = 0, r G Z. Here uo is the random variable which takes a 
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value in N. u{r, t) = if t — r is even or if |r| > t, and 

w(i,j)=uj(i-j,i+j-l), (2.2) 

G are geometric random variables. In our model, all w(i,jys are not independent 
as in 1.19 , but a symmetry condition w(i,j) = w(j,i) is imposed. The parameters of the 
independent geometric random variables are given by q (resp. 'jy/q) for off-diagonal (resp. 
diagonal) points, 

F[w(t,j) = k] = (l-q)q\ l<j<i, (2.3) 
F[w(i,i) — k) — (1 -l^q){l^q)\ 1 < i, (2.4) 

for fceN. Introduction oiw(i,j) above is almost unnecessary for the following discussions, 
but useful to see the connection to the problems of last passage percolation and ASEP 

In the following, we devote our most discussions to the two special choices of 7 where 
the detailed analysis of the model is possible. One choice is 7 = 1, which is known to 
correspond to a critical point of the model j2H|. The other is 7 = 0, which corresponds 
to the model without the external source at the origin. In the following, the 7 = 1 (resp. 
7 = 0) case will sometimes be referred to as the orthogonal (resp. symplectic) case because 
it will turn out to be related to the orthogonal (resp. symplectic) to unitary transition of 
symmetries in random matrix theory. The analysis of the two cases are quite analogous, 
so that we mainly consider the 7 = 1 case in the sequel and summarize the results for the 
7 = case in sectional 

We would like to know the statistical properties of h(r, t). As already mentioned in the 
introduction, the fluctuation of the height at a single point is already known. It is given 
by the GOE Tracy- Widom distribution at the origin j2H] and by the GUE Tracy- Widom 
distribution at other points [7j. In this article, we are interested in the crossover between 
these two. Following [T%lll4[ll9j . we introduce the multi-layer version of the model which is 
defined as follows. In the original single-layer discrete PNG model, when two steps meet, 
the higher step takes over the lower one. The information of the lower step is lost. To keep 
this information, let us suppose that there was a second layer below the original layer from 
the outset. The height of the second layer remains —1 until a collision of steps happens 
at the first layer, at which time the nucleation for the second layer occurs just below the 
collision. The height of the nucleation is taken to be the height of the overlapped region 
of the higher and the lower steps which collide at the first layer. See Fig. 1. The steps 
in the second layer grows laterally in both directions with unit speed as in the first layer. 
The nucleations at the second layer are only due to collisions of steps at the first layer; 
there is no additional probabilistic nucleations. We denote the height of the second layer at 
position r and at time t as hi(r, t). One can further suppose that there were infinitely many 
layers below the first layer, equally spaced at time and successively construct hi(r,t) for 
% = 2, 3, • • • . We also set h(r, t) = ho(r, t). This is the multi-layer PNG model. An example 
of the whole height configurations is given in Fig. 2. 
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In the following, we mainly consider an even time, which is fixed to t = M = 2N 
until the asymptotics is considered. There are infinitely many lines, but for time t = M 
the number of lines which are not straight is at most N. Moreover, there is a restriction, 
hi(M,M) = = 0, 1, • ■ ■ , N — 1). For each jump with length k of all layers, we 

associate a weight (*Jq) k - Since the whole configuration of the layers keep all information 
of nucleations, the product of all weights of all jumps in the multi-layer PNG at t — M is 
given by 

N 

II q^HiVqT^. (2.5) 

i+j<2N+l,i>j i=l 

From now on, we change the way of looking at multi-layer heights. The space coordinate 
r in the original setting will be interpreted as the time coordinate. The height coordinate 
will be interpreted as the space coordinate and is represented as x. We set hi(r,t) = x r i+1 
for % = 0, 1, • • • , N — l and r = 0, 1, • • • , M, and regard the height lines as paths of particles. 
Let us denote the collection of the positions of particles at time 
where x\ > ■ ■ ■ > x r N , and a collection of them for all times as x = (x°, • ■ ■ , x M ). 

In what follows we use a = ^Jq. Suppose that, when the the particle picture of heights 
are employed, the x and r axes are pointing the right and the up respectively. Then at 
odd times, particles move to the right. The transition weight of a particle from x at time 
2u to y at time 2m + 1 is given by 

^,^v) = \ ( y a)a '"' v ~ x ' (2.6) 

[0, y < x. 

At even times, they move to the left. The transition weight for each particle is 

a i \ J ' y> 
I (1 — at)a x y , y < 

Then, according to the Karlin-McGregor theorem [2111110] > the weight of non- intersecting 
paths of n particles which ends at xf 1 = 1 — i is given by 

M-l 

w nM {x) = n det(0 r , r+1 (x[,^ +1 ))^ =1 , (2.8) 

r=0 

with the restriction xf 1 — 1 — The weight is written in a determinantal form, and hence 
the corresponding process is called a determinantal process. Notice that the weight in ()2.8|) 
with n = N is slightly different from the weight of the multi-layer PNG model because 
some configurations contained in ()2.8|) are not allowed in the multi-layer PNG model where 
there are infinitely many particles. On the other hand, if we take the limit n — > oo from 
the beginning, there would appear another difficulty due to the infinity of the number of 
particles. Thus, for simplicity of the treatment, we approximate the weight of the multi- 
layer PNG model by (|2.8J) for the moment and take the limit n — >• oo afterwards when it 
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(2.7) 



becomes easy. The simple determinantal expression in (|2.8|) allows us to perform a detailed 
analysis of the model. If we took diagonal parameters to be more general as in (|2.4jl . the 
expression would become more complicated. 

For later use, we introduce the generating function of <fi r ,r+i, 

Mz) =Y,<t>2u,2u+i{k)z k = (2.9) 
z — ' 1 — az 

fcez 

Let us denote the Fourier coefficient of a function a(z) as a(k) like 



][>(*)**. (2-11) 



For instance one sees 



In addition, for a product, 



the Fourier coefficient, 



a{z) = 



!>2«,2«+i (x,y) = f 2u (y - x), (2.12) 
^2u-i,2u(x,y) = hu-i{y - x). (2.13) 

I '2 

fn,r 3 (z) = l[f l {z), (2.14) 

l=ri 



</>n, n {x, V) = fn,n(y ~ x )> (2-15) 
represents the transition weight of a particle between (r 1; x) and (r 2 ,y). Notice that when 
r i = r 2 — r one has 

4>r t r(x,y) = 5 xy . (2.16) 

Before going to the next section, a remark is in order. Although we have employed the 
situation where the geometric random variables are given by a single parameter q, up to 
some point, it is easy to generalize our discussions to the case where 

F[w{i,j) = k] = (l-a i a j ){a i a j ) k , l<j<i, (2.17) 
F[w(i, i) = k] = (1 - di)a k , 1 < i, (2.18) 

with < a* < 1 for all i. Then (f271.(|2~T7]il are replaced by 

/2u(«) = i , (2-19) 

1 — CLN+u+l z 

f2u- 1 (z) = 1 , (2.20) 

1 — ClN-u+l/ Z 

with the transition weight of a particle again given by (|2.12j) . (|2.1Hj) . It should be noted that 
our free parameters are only a/s; for the case in |18j . there are two sets of free parameters 
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3 Determinantal Process and Kernel 

Let us consider the weight of non-intersecting paths given by 

M-l 

WnM*) = II det(0 r , r+1 (x[, )V; J ,. (3.1) 

r=0 

where xf (i = 1,2, ••• ,n) is fixed. As for there is no restriction. This is the main 
difference between our model and the model in ^J] i n which x® is also fixed. The matrix 
element r . jr+ i(x[, a^ +1 ) is the transition weight of a particle between (r, x\) and (r+1, a^ +1 )- 
Our main focus is on the special case in (|2.fij) . (|2.7|) . but the results of this section do not 
depend on a specific choice of 0. For instance, when particles can hop only to nearest 
neighbor sites for each time, the determinantal process is nothing but the vicious walks 
studied in [HUES EH]- Hence our discussions below can immediately be applied to the 
problems as well. 

The partition function is defined as 

Z n,M = JwE^^' ( 3 ' 2 ) 

x 

The probability of the non-intersecting paths reads 



1 



Pn,M( X ) = f „y\ M r7 W nM X )- ( 3 - 3 ) 



Let g(r, x) (r — 1, • • • , M) be some function. In this section, we show 
Proposition 3.1. 

M-l n 

E II + 9{r,x^)p nM {x) = y/det(l + K l9 ), (3.4) 

x r=0 j=l 

where the determinant on the right hand side is the Fredholm determinant, 

oo j M M 

dfit(i+%) = EjfciEE E-EEE g(n,x 1 )--. g (r k ,x k ) 

k=0 ' n=l xi ji=l,2 r k =l x k j k =l,2 

det {K x (ri,xi\ r v , x v )j l ,j l ,)l v=1 , (3.5) 
and K\ is the 2x2 matrix kernel, 



Ki(r 1 ,x 1 ;r 2 ,x 2 



S 1 (r 1 ,x 1 ;r 2 ,x 2 ) D x {r u x x ; r 2 , x 2 ) 
I 1 (r 1 ,x 1 ;r 2 ,x 2 ) S , i(r 2 ,x 2 ;ri,Xi) 



(3.6) 
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with the matrix elements being 

S 1 (r 1 ,x 1 ;r 2 ,x 2 ) = S 1 (r 1 ,x 1 ;r 2 ,x 2 ) - n , r2 (xi, x 2 ), (3.7) 

n 

Si(ri,xi;r 2 ,x 2 ) / <t> ri ,M{xi, x f)( A i 1 )i,j ( ^i( r 2, x i\ M , x f), ( 3 -8) 
h(r 1 ,xi;r 2 ,x 2 ) = Ii{rx,xi;r 2 ,x 2 ) - G l (ri,x 1 ;r 2 ,x 2 ), (3.9) 

n 

h(n, x x - r 2 , x 2 ) = - x u M, x^){A-\p x {r 2 , x 2 ; M, xf ), (3.10) 

n 

D 1 (r 1 ,x 1 ;r 2 ,x 2 ) = E n ,Af(zi, ^fX^ki^M^, x* f ). (3.11) 

ifere 

(Ai)ij = E s#n(y 2 - yi)(f)o,M(yi, xf 1 )(f) 0i M(y2, xf), (3.12) 

2/1,2/2 

Gi(ri,xi;r 2 ,x 2 ) = E s#n(y 2 - yi)(f)o, n {Vh Zi)0o,ra(j/2, x 2 ). (3.13) 

J/1,2/2 

T/ie expression for general Ki(ri,Xi;rj,Xj) (i,j = 1,2, ••• ,M) has the same form as 

Remark. The subscript 1 refers to the fact that this case is related to the orthogonal- 
unitary transition in random matrix theory. If we take g(r, x) = —x.J r i x ) ( r = 1> 2, • ■ • , M), 
where \j r is the characteristic function of J r , ()3.4)1 gives the probability that there is no 
particle on J\ x ■ ■ ■ x J M . It should also be remarked that at this stage the infinite particles 
limit is taken easily. One only replaces the summation, YLij=n m (|3.8|> . ()3. 10|l . ()3. 1 If) with 

Proof. We derive (|3.4f) by generalizing the methods of ^§1 and [31]. Let us define 

Af-l n 

Z n , M [g] = E II II^ 1 + ^)H,m(x). (^) 

V ' x r=0 j=l 

Remark Z* nj jvf[0] = Z Uj m so that 

n 

\~f I I -J- n(r rr r \\r> . A ' v\ — . 

^n,Af[0] 



M-l n 7 \ 1 

e n n^+^^^w^) = (3.i5) 



x r=0 j'=l 

By repeated use of the Heine identity, 



i E det(^(x,))^ =1 det(^(x J ))^ =1 = det ( E &0*>iPO ) , (3-16) 



n! 



s 



one has 

ZnM= E det ( E ^ + 9(0,x i )) ( j>o, 1 (^,X 1 )(l + g(l,X 1 ))--- 

x°>—>x% Xi,--,Xm-1 

(X M _ 2 , X M -i)(l + tf(M - 1, X M -i))0m-i,m(X M -i, xf ) ) . ( 3 -IT) 
Now, using another identity 



' ' i,j=i 



E det (^(^))L=i) =det ^sgnfe-^^y , (3.18) 
one finds 

Z n ,M\g] 2 = det^E sgn(?/2 -yi) E ( X +5 , (°'?/i))^o,i(l/i,Xi) • • • 4>m~i,m( X m-i, xf) 

2/1,2/2 Xi,---,X A /-i 
E ^ + fi '^ ' ?/2))0O,l(2/2, ■ • • 0m-i,m(^m-i, ^f) 



»J=1 



= det(E s g n (?/2 - J/i) 

2/1.2/2 

E 00,0(2/!, *o) (1 + ^(0, X ))MX , X x ) ■ ■ ■ 0m-i,m(^m-i, af ) 

-Xq,"' ,-XiW-l 

E Mjfe, X )(l + <?(0, X ))MX , X x ) ■ ■ ■ 0m-i,m(^m-i, xf ) 



(3.19) 

For the second equality we used (f)ofi(x, y) = S x>y . We divide the element of the matrix 
in the above determinant into four parts with or without factors of the form g(r,X r ) or 
g(f, X r ). We have 

Z n , M [gf = det((A 1 ) 4J + (4 1} )m + ( A i ] h + ( A ?)i,j) n , (3-20) 
where (Ai)i t j is already defined in ()3.12|) and 

M 

(^ki = E sgn ^ 2 ~ E E E 0o,n(yi,^i) 

2/1,2/2 1=1 0<n<-<ri<M Xi,- ,X t 

1-1 

l[g(r s ,X s )<p rs , ra+1 (X s ,X s+1 ) ■ giruX^uiXuxf) ■ <f> 0M {V2,xf), (3.21) 



M 



(4 2) )m = Y sgn ^ 2 ~ foMvi'Xi 1 ) Y Y Y ^Av^^i) 

J/i,J/2 1=1 0<ri<-<f t <M x 1 ,- ,X t 

l-l 

]Jg(r s ,X a )<l> f ., fa+1 (X a ,X a+1 ) • gifuX^MiXu^f), (3.22) 



s=l 

M 



(4 )<j = Y sgn ^ 2 ~ ^) Y Y ^ Xi ) 

2/1,2/2 «=1 O^K-^rKMX!,-^ 

r s ,r s+ i(X s , X s+ i) ■ g(ri, Xi)(j) rh M(Xi, x i ) 

8 = 1 

M 

Y Y Y ^0,^(2/2,^1) 

m=l 0<fi<-<f m <M Xl ,- ,X m 
m—1 

JJ #(^,^)0r s ,f s+1 (^ S ,^ S +l) • g(fm,X m )(f)f miM (X m ,xf ). (3.23) 
s=l 

Let us define 

¥»(ri,xi;r 2 ,x 2 ) = S^i, xi)0 rijra (xi, x 2 ), (3.24) 
^(ri,xi;r 2 ,x 2 ) = 0n,r- 2 (^i, x 2 )g(r 2 , x 2 ), (3.25) 

and 

(f>i(r,x) = (p r>M (x,xf), (3.26) 
(Gi) i (r,x) = Gi(r,x;M,xf f ), (3.27) 

^ifrs) = Y Y Y <Pr,r 2 (x,X 2 )Y[ip(r s ,X s ;r s+1 ,X s+1 ) ■ ip(r h Xf,M,xf*). 

1=1 r 2> - ,n X 2> - ,Xj s=2 

(3.28) 

Then (A^)jj can be rewritten as 

M M-l M-l 

(A?) id = YY---Y Y G 1 (r 1 ,X 1 ;M,xf)g(r 1 ,X 1 ). ( f )rur2 (X 1 ,X 2 ) 
1=1 n=o n=i 
2-1 



JJ tp(r s , X s ; r s+1 ,X s+1 ) ■ tp(r h Xf, M, xf 

s=2 

= Y9(r,x)i; i (r,x)G 1 {r,x;M,xf) 

r,x 

= YaA-{Gi) r (3.29) 
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In the last expression, the dependence on r, x is omitted for notational simplicity. Noticing 
the antisymmetry of G\, 

G 1 (r 1 ,x 1 ;r 2 ,x 2 ) = -G 1 (r 2 ,x 2 ;r 1 ,x 1 ), (3.30) 
we can also rewrite (Ai)\j and (Ai)\j as 

(A? } ) hJ = -J2g^-(Gih (3.31) 

r,x 

= - ^^gtpj ■ Gi(gipi), (3.32) 



r.x 



where we used the notation, 

(Gif)(r, x) = G i( r ' x ' r i> x i)/( r i> x i)- ( 3 - 33 ) 



Hence we find 



ZnMQ? = det ^(v4i)jj + YM&i ■ - 9^ ■ {Gi)i - ^ ■ G 1 {gijj i ))\ . (3.34) 

V r > x ' i,j=i 

In view of (J3.15)) . we divide this by 

Z nM [0] 2 = det A v (3.35) 
In the determinant this corresponds to multiplying by A^ 1 , say from the left. Let rji = 

ZMi 1 kA^i = EMi 1 ki( G Jr Then 

(fjfjoj) = det + • • ( 3 - 36 ) 



One notices that, if one introduces 



B(i;r,x) = (-g^i- gGi(gr]i),gr]i), C(r, x; i) = ( \ 1 , (3.37) 



this is equal to 

det 1 8 id + J2B(i;r,x)C(r,x;j)) = det(l + BC). (3.38) 

V r, a / iJ=1 

Now we use a simple fact 

det(l + 5C) = det(l + OB). (3.39) 



11 



The determinant on the left hand side is the determinant of a matrix of finite rank, whereas 
the one on the right hand side is the Fredholm determinant. We see 



det(l + CB) = det 



i - E* A ® 9®i -T.i^i® gGi(gvi) E* V>< ® gvi 
- Ei(Gi)< ® - Ei( G 0i ® ^Gi(^) 1 + £i(c?i)< <g> 



where we used the notation a®b for the operator with kernel of the form a(r\, X\)b(r2, x 2 ) ■ 
The matrix in the determinant may be rewritten as a product of two matrices, leading to 



det 



1 - Et i>i ® Ei V>i ® fl 1 ^ 



1 0' 



(3.40) 



- Ei A ® Ei V 1 * ® 0*71 

-Ei(Gi)i®^i-Gi 5 J2i(Gi)i®9Vi 



M 



The determinant of the right matrix is one, so that this equals 
det ^1 + 

Let us remember the definition of ipi in (J3.28)) and notice 

= (1 - vY 1 * 4>» 

g^i = * = (1 " ¥>) _1 * (9<t>i), 

1=1 

where 

(/1 * f2)(ri,x 1 ;r 2 ,x 2 ) = ^2fi(r 1 ,x 1 - 1 r,x)f 2 (r,x- 1 r 2 ,x 2 ) 

r,x 

The determinant we are considering is now 

-Ei&®0*< Ei^®^(^rV)i 



Af 



0+ 


"(i-yT 1 o" 






1 





Multiplying from the left by 



and from the right by 



1 = det 



1 = det 



1-0 
1 

1 

1 - V 



(3.41) 



(3.42) 
(3.43) 

(3.44) 



1 

(1-vr 1 

(3.45) 



(3.46) 
(3.47) 
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one gets 



det 



l-Cp 

o i - V 



+ 

det I 1 + 



5 



i — i \ 



E*.,(Gi) 4 ® (VM^ - Gt E l AG 1 ) l ® (AI 



(3.48) 

Recalling the definitions of <pi and (Gi), in (|3.26j) and (|3.27|) . we see that (|3.4j) holds. □ 



4 Scaling Limit 

In this section we consider the scaling limit, where the universal properties of the model 
are expected to appear. 



4.1 Generating Functions 



To study the asymptotics of the kernel, we compute the generating functions of Si, Ii,D\. 

First let us recall some basic facts about a Toeplitz matrix jSH]. A Toeplitz matrix is a 
matrix A of the form, 

CLq CL—i CL—2 ' ' 



A 



a 2 cl\ 



«o 



It is useful to define a function, 



a k z 



(4.1) 



(4.2) 



which is called the symbol of the Toeplitz matrix A. Conversely, for a given function 
a(z), set dk = a(k), the Fourier coefficient of the symbol a(z). Then one can define the 
corresponding Toeplitz matrix (|4.1j) . which will be denoted by T(a). The inverse of A 
can be given in a compact fashion, if the symbol a(z) has winding number zero and is 
Wiener-Hopf factorizable as 

a(z) = a + (z)a_(z) (4-3) 
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with 

fc=0 

oo 

Y,a-{k)z~ k . (4.5) 



oo 

a+(z) 

fc=0 



a-{z) 

k=0 



In fact for a factorization ()4.3|) one has 

A = T(a) = r(a_)T(o+). (4.6) 
In addition, the inverse of the matrix T(a±) is simply given by 



(4.7) 



Hence the inverse of the matrix A is 



A- 1 = T ( — J T f — ) . (4.8) 



a + J \a 

To obtain the generating functions of Si,I\,Di, we need the inverse of the matrix A\ 
(|3.12|) . As remarked below ()3.13|1 . for the multi- layer PNG model, the matrix A\ should 
be taken to be an infinite dimensional one, which turns out to be a Toeplitz matrix. If we 

set 

{A 1 ) ij = a{i-j), (4.9) 

the symbol is computed as 

a x {z) = y^ j a(k)z k 
fcez 

= kM(l-k-y 1 )z- 1+k+ ^-sgn(y 2 -y 1 )zy*-yi-f , M (l-y 2 )z 1 -y> 
k,yi,y2 

= fo,M (-) Sl (z)f 0M (z), (4.10) 



where 

Sl (z) = Y j sgn(h)z k = l -±^. (4.11) 
^— ' 1 — z 

In the second equality of (|4.1()j) . we use (|2.15j) and (|3.12j) . The winding number of 01(2;) is 
not zero and hence the formula (|4.8j) does not hold as it is in this case. But the difficulty 
can be overcome fairly easily. Let us first notice that s\(z) can be written in two ways as 

Sl (z) = S U(z) S U(z) = sfl(z)sfl(z), (4.12) 
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with 

45w = ^, 4 1 2w = i, (4.i3) 

s «(z) = 1, sfl(z) = -l±ig. (4.14) 
Correspondingly, if we define 

af + { Z ) = / ,M;- W(«)4W» ( 4 - 15 ) 

afl(z) = /0,M;+ (4.16) 

for i = 1,2, we have 

ai(z) = ag(^)a^(z) = afl{z)afl{z). (4.17) 
Then it is confirmed that the inverse A^ 1 is given by 

v -H T (i) T fe) +r (i) r (i)^ (418) 

which is the substitute for ([4.8)1 . 

Now one can compute the generating function, 



fe=i I iez jez 

+ £ -r +fe ( v^S) c< - k) y: (iaS) (* - 3) 



1 



fc-1 



+ 



_ __Z1 1 

Z\ - Z 2 fo,M;-{Zl)fo,M;+(-^)fo,M;+( Z 2)fo,M;-(~^) 

x n4!(i)W) + 'S(^W)}' (419) 

which is valid for |zi| > j^. In the first equality we use a fact that (l/ Q i+) ( k ) = for 
k < and (l/ a i-) (A;) = for k > (z = 1, 2). 
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We also introduce the generating function of the function G\ in ()3.13|) . Note that as a 
function of Xi,x 2 , G\ only depends on the difference Xi — x 2 , so that we can set 

G?i(ri, xi; r 2 , x 2 ) = (G ! i) ri>ra (a: 1 - x 2 ). (4.20) 

Then one has 



(Gi) ri ,r 2 (z) = ^2{G 1 ) rit r 2 (k)z k = /o,n(«)si (-) /o,ra (-) = -/^(^lW/o, 



1 

z y 
(4.21) 

With these preparations, calculations of the generating functions of Si,Ii,Di are not 
difficult. We have 

"2 



S 1 (ri,z 1 ;r 2 ,z 2 ) = Si(r 1 ,xi;r 2 ,x 2 )z* 1 z 2 

"EE ^1^(1 - « - Xl)^ 14 * 4 " 1 ■ 4" i (^r 1 )i J -4~ 1 • (Gl)r 2 ,M (X 2 - 1 + j) Z 2 
xi,xa i,j 



zx fn,M--(^)fo,r2',+(^)fo,M;-( z 2) 1 / S l+(^) s i+(^) 



(1- 0^)^(1 -o/za)^ 21-^1 (1 + Zl)(^2-1) 

ii(ri,zi;r 2 ,z 2 ) = ^ ^( ri > ^ r2 > x 2)zi 1 z 2 X2 

- —(Gi) ruM {zi) ■ y^ z l ' 0*4-1 l )i,j z 2 1 ' (Gl)r 2 ,M ( — 
i,j 

Zx /o,ri;+(^l)/o,M;-(^-)/o,r 2 ;+(^)/o,M;-(2 ; 2) 
z l~ z 2 fri,MA Z ^)fr%,MAh) 



1 ^ ,(1) I 1 ^ „(1) ( l \ > J 2 ) ( 1 ^ J 2 ) ' 1 



(l-o) 2 ( ni+U2 )(l-a/2 1 ) 7V - Ul (l-o2 2 ) A '- M2 21 (lz 2 + l 11 + 



~i 



(1- 0^1)^1(1 -a/z2) N+ua z 1 -z 2 \2z 2 -l 21- Z! 



(4.23) 
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h{ri,zx;r 2 ,z 2 ) = 



^2 D(r 1 ,x 1 ;r 2 ,x 2 )z Xl z 2 x 

— fri,M ' 1 )i,jz{ ■ f\ 



Z\ 



r 2 ,M\ Z V 
fr 1 ,M;-('k)fr2,M--i Z ^) 



" Zl~Z 2 /o,M ; -(^)/o iri;+ (i)/0 >r2 ; + (^)/, 2> M ; -(i) 2 \ *<2(£)*£(£) 8${±)s?l{±) 

(l-a/z 1 ) N+u ^l-az 2 ) N+u * z x \\ z x -\ 1 1 - z 2 \ 

(1 - a) 2 ('"i+«2)(l - azi) N ^*(\ - a/z 2 ) N ~ u * z 1 - z 2 [2 zt + 1 2 1 + z 2 J ' ^' ^ 

In the last expressions of these, we have set r, = 2itj for i = 1,2. 
4.2 Bulk 

We are now in a position to study the asymptotic behaviors of the model. The thermody- 
namic limit shape is already known to be 



h(r = 2P N, t = 2N)/N ~ [a + , (4.25) 



2a 
1 — a 

where < (3$ < 1 is fixed [37J. As for the fluctuation around this limit shape, the KPZ 
theory tells us that the correlation survives for 0(N 1 ^ 3 ) along the height direction and 
0(N 2 ^ 3 ) along the r direction. 

Let us define the scaled height variable as 

„ , a , Hr = 2(3 N + 2c(Po)Nlr, t = 2N) - q(fl, + ^)N 

Hn{t,(3 ) = , 4-26 

where 

a ^ = YT^2 { a + V^P) , (4-27) 
i 

d W = 7, ,w! ^^ V^l+^g + a V / l^) i (v / l^ + + (4.28) 

(1 — a 2 )(l — p 2 )6 

c{(3) = a"3(l - /3 2 )t(v/l + £ + ay/1 -(3)^{^Jl-(3 + ay/l +P)*. (4.29) 

Remark that a((3) just comes from the limit shape ()4.25|) whereas d((3) and c({3) are taken 
so that the results become simple. In this subsection, we show that the fluctuation of the 
model in the bulk is described by the Airy process f2JEI]. Namely, we show 

Proposition 4.1. 

lim F[H N (n, A)) < si, • - • , H N (r m , fa) < s m ] = det(l + K 2 Q). (4.30) 

N^oo 



17 



The determinant on the right hand side is the Fredholm determinant, 

oo m » m „ 

k=o ' n=l -r fc =i 

det(/C 2 (r^/;^,^))t'=i' ( 4 - 31 ) 

where G(Tj,£) = —X(sj,oo)(%) (j = 1,2, ••• , m) and ]C 2 is a scalar kernel, for which a 
representative element is given by 

/C 2 (n, 6; 72, 6) = Wn, 6; r 2 , 6) - ^ (6, 6), (4.32) 



/•oo 

^ 2 (Ti,6;-r 2 ,e 2 )- / rfAe^ Tl - T2 )Ai(ei + A)Ai(6 + A), 
Jo 



(4.33) 



■ 5 T\ ^ T2 j 

<,, ^-^'--' ) <I _n o rfAe-^--)Ai(e 1 + A)Ai(6 + A), n < r 2 . (434) 



Remark. The kernel ()4.32j) is called the extended Airy kernel and had already appeared 
in the context of Dyson's Brownian motion model jSHUSH]- We also remark that some 
properties of the Airy process have been discussed in jTUHniliniHIl • F° r the special case 
of a single time, we have 

lim F[H n (t,/3o)<£] = F 2 (Z), (4.35) 

N— +00 

where F 2 (£ ) is the GUE Tracy- Widom distribution [20] • In the language of the last passage 
percolation, (|4.35|) means that the fluctuation of the last passage time at an off-diagonal 
point is described by the GUE Tracy- Widom distribution. For the nonsymmetric case, the 
statement has been already established in But for our symmetric case, this has not 
seem to be shown although it was already expected based on a physical argument in 
Remember also that the fluctuation at an off-diagonal point has not been handled by the 
techniques using the Riemann-Hilbert method j25| . 

Proof. Let us derive (j4.3Up . We first remark that ()4.33|) has a double contour integral 
representation, 

Iff e i(iw-i+i&w 2+^(wf+w%) 
Ktffa, fi; 72: &) = --— / dw X / dw 2 —7 ; r, 

when 771 + r] 2 + Ti — t 2 > 0, 771, r] 2 > 0, and that ()4.34|) can be computed as 

$ njT2 (^ 2 ) = 1 e ^-^mH^-nr (4.36) 

A /47r(r 2 - n) 
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when r 2 > T\. These are already given in [Hi] . 

Using the generating functions in the previous subsection, we get the double contour 
integral formula for Si,I\,Dx, 

S 1 (r 1 ,x 1 ;r 2 ,x 2 ) = 1 I ~^ff dz 2 z 2 2 ^ 1 {r 1 , z x \ r 2 , z 2 ), (4.37) 



(2™) 2 Jc Rl zT +1 J c 



R 2 



I 1 (r 1 ,xi,r 2 ,x 2 ) = ,„ / -zSr / dz 2 z 2 2 r li(rx,zi;r 2 ,z 2 ), (1.38) 



(2^) 2 Jc Rl zT +1 Jc R2 
Di(r 1 ,x 1 \r 2 ,x 2 ) = , * 2 / -ffrr / fe2 2 ~ lro i( r i>^i; r 2,^)- (4.39) 

Here C# denotes a contour enclosing the origin anticlockwise with radius R. We have taken 
Ri > R 2 to agree with the computation in (j4.19j) . For Ii, exchanging the radius of the 
contours of z\ , z 2 gives 

Ji(ri,xi; r 2 ,x 2 ) 

%l dz 2 zr X Uru^r 2 ,z 2 )-^\ dz 2 zr^\G^J^\ 
m ) Jc R „z 1 J Cri 11X1 Jc R2 \ z 2/ 

dz 2 z 2 2 ~ 1 I 1 (r 1 , z x\ r 2, z 2 ) + Gt(ri, xi] r 2 , x 2 ). (4.40) 



(2™) 2 Jc R2 zT +1 Jc Rl 
1 f dz\ 



( 2 ™) 2 Jc R2 z? +1 J CRl 

In the first equality the second term on the right hand side appears from the contribution 
of a pole at Z\ — z 2 . Hence one gets the double contour integral representation for ii, 

1 f dz f 

h{ri,x 1 ;r 2 ,x 2 ) = / -— ^ / dz^^hin, zi,r 2 , z 2 ). (4.41) 

{2my j Cr2 Zl 1+ j Cri 

Let us first consider the asymptotics of Si. Since we can follow the discussions in |l9j 
for a large part, we omit details and state the main steps. If we write 



we have 



g^(z) = (1 + (3) \og{z -a) - (1 - /?) log(l - az) - + /3) log z, (4.42) 



S 1 (r 1 = 2u u xi\r 2 = 2u 2 , x 2 ) 

( 1 _\2(« 2 -«l) r Av r rlv -,352-^(^2-1) 



(2vrz)2 J Cri zi J CR2 z 2 ^-"(mi-D 
x-3— Zl ~^ 2 J , (4.43) 

^1-^2 I (l + Zl)(z 2 -l)J ' 1 ^ 

where (3i = Ui/N, j3 2 = —u 2 /N. fii, \i 2 are arbitrary constants at this stage. We would like 
to apply the saddle point method to this integral. In general for each value of fi, there are 
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two critical points zf = p(fi, j3) ± ^/p(p,, f3) 2 — j3) with 

qfaff) = £±l. (4.45) 
\i p 

Since we are considering the scaling in ()4.26|) . we will set 

x i = Na(0 i ) + d(J3 o )N*Z i , (4.46) 
with £i the scaled space variable. Accordingly we set /ij = /i c (A) with 

fji c (j3) = a(/3) + 1 = (l + « 2 + 2a v / l-/5 2 ) • (4.47) 

For this special value of /i, the two critical points, £jr, merge to the double critical point 
PM, where g'^^pdP)) = g'^pMP)) = 0, given by 



VI - p + ayl + P 

The main contribution to the integral in ()4.43|) comes from around the double critical 
points z\ ~ Pc(Pi), Z2 ~ l/pcifti)- The paths of integration may be deformed to 

Zl = PciPl) (l ~ M aLrU^ ) > ( 4 - 49 ) 



rf(A)iVV3 

- = p c (/3 2 ) (l - —l—- w 2 ) , (4.50) 



z 2 n ™ \ d{p Q )N 1,Z 

where wi = Q + irji with — oo < Q < oo (i = 1, 2) and 771, 772 > are fixed. For a fixed /3q, 
we take 

r, = 2m = 2N U + C ^0} , (4.51) 
where r, is the scaled time variable. To leading order we have 

Pc(Pi)=Pc(Po)+p' c (Pom-Po)=Pc(Po) (l+ » ( 4 - 52 ) 



Pc(Jh) = Pc(-A)) + Pc(-/3o) (& + A)) = Pc(-A)) ( 1 - 37^777? ) > ( 4 - 53 ) 



rf(/5 )iW3 

where 

lPoJ rfM)d{fa) p' c {-W%) ' 

is used. 
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Expanding around Z\ ~ Pc(Pi) to order 0(1), one finds 



^ c (/3i)A(Pc(/5i)) + ^? (4.55) 



after some computation. Similarly one gets 



Ar ^ c (/3 2 ),/3 2 ( 1 /^) ~ Ng M)> p 2 ip c ((3 2 )) + -«4 (4.56) 

Let us denote 

Ao(^) = 9^(p)AVc{l3)) 

= 2/31og(l - a 2 ) + 1(1 + 0) log(l + /3) - 1(1 - 0) log(l - (3) 

+ T-^—ziy/l-P + a^/TTp) 2 log(Vl-^ + a^TT^) 
1 — a z 

- — !-j( v/TT^ + av^) 2 \og(V^T0 + a^/T^P). (4.57) 
1 — ar 

Clearly, this is an odd function. Therefore if we expand as 

X(J3) ~ A (AO + A!(A))09 - A)) + A 2 (/W - /^o) 2 + A 3 (/5 )(/3 - A)) 3 + • • • (4.58) 
we have Aj(/3o) = ( — l) i+1 Ai( — /?o) for i = 0, 1, 2, • • • . Expanding to order 0(1), we have 

Ng^MPi)) ~ iVA (/?o) + A 1 (/? )c(/?o)iV 2/3 r 1 

+ A 2 (/3 )c(/5 ) 2 iV 1 / 3 r 1 2 + A 3 (/5o)c(/3 ) 3 r 1 3 , (4.59) 
Ng^m^M) ~ -WAo(ft) - A 1 (/5 )c(/5 )iV 2 / 3 r 2 

- A 2 (/3 )c(/3 ) 2 iV 1 /3 r 2 _ A 3 (/5o)c(/3 ) 3 r 3 . (4.60) 

Combining dUMD , dHHOD with (1432J) . (14331) gives 

^ = ^( 1 + 5(^V3( T '- imi 0- < 461 > 



so that 



x 2 +V(l-^ c (/3 2 )) 

fp(p\y&-ti)d(f3 )N 1 / 3 e &T2-ZiTi+i€iwiMbm_ (4.63) 



si+iV(l-/ic09O) 

i 
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In addition, to leading order we have 



* r i+ /I ., (4.64) 



Zl-22l (1 + ^l)(^2 — 1) J T"2 - Ti + + W 2 ) 

to which the second term on the left hand side does not contribute. 
Finally one obtains 

§1 ~ (1 -^^-"^Pc^o))^^^^ 3 ^^ 



e 



Ai(/3o)c(/5b)Ar 2 / 3 (Ti-T2)+A 2 (A))c(A,) 2 ^ 

1 ; dw x [ dw 2 ( , J e «i«i-H6««»+i(«'?-H»S). (4.65) 



4 ^ 2 Jlavn^ Jjmw 2 = m V r 2 - ri + z(wi + w 2 ) 
We also need the asymptotics of 

(j) ri:r2 (xi,x 2 ) 
(1 - af (u2 - u ^ f dz 



2m 



f dz 

/ —z X2 - Xl [(l-az)(l-a/z)] U2 ~"< 
JCi z 



_ (1 a) 2{ - U2 Ul > I ^^-^(/^-lj-^+jv^ {AM) 
2iri J Cl z 

Let us set 

From 1)4.52)1 . this can be rewritten as 

Expansion around z ~ p c (Pi) as in (|4.55|1 gives 

Ng^faWxW ~ ^08 1 ),/Ji(Pc(A)) " 3^1 - i*) 3 . (4.69) 

Similarly 

;~^)( 1+ 5(A^ (T2 - i<r) ) (470) 



leads to 



-V//,,,...^...^ ( -) ~ N gt , cmA ( Pc ((3 2 )) + i(r 2 - «t) 3 . (4.71) 



One also finds 

^ 2 -AT( Atc (,8 2 )-i)- a; i+iV(Mc(A)-i) ^ (p c (^))^)^ 1/8 (6-&) e M&-&). 
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Hence one gets 

^^(xx.xa) ~ (1 - a) 2 ^-^)(p c (/3 )) fe -^ )ci(/3o)iVl/3 



e A 1 (A,) C (A,)W 2 /3(ri~r 2 )+A 2 (/3o)c(/3o) 2 AfV3( T 2_ T |) + A3( A) ) c ( / 3 )3(r3-r|)-^ + ^| 

— / dae^"^-^-^-^ 2 . (4.73) 

For 7i, one has 

7i(ri = 2u 1 ,x 1 ;r 2 = 2u 2 ,x 2 ) 



f dzi f 



(2 m y J CR2 Zl J Cri Z2 ^i-^o-u 



e N(9 tlc(M ,f3 1 (zi)+g tlc{l32h g 2 (l/z2)) 



X z l -z 2 ^2z 2 -l 21-^}' (4 ' 74) 

where $ = —Ui/N. The critical points of 21,22 are the again given by 21 = p c (Pi),z 2 = 
l/p c (P 2 ). To leading order, the paths are 



^=^\ l -im^ {Tl+iWl) )' (475) 

ir M - M ( l -4m^ iT2+ ' W2) )- (476) 



so that we have 



Jx ~ (1 - a) 2(Ml+U2) 7 . |4t7^ 7 .^J Pe(/3o)) (gl+6)d(Al)JVl/3 - 



(1 + Pc (-A)))(i - p c (-/?o)) 3 v " v " 77 d(A,) 3 iv 

/ dtu! / du> 2 (ri - r 2 + - w 2 )) e <&»i^+i(»M). (4.77) 

</lm«)i=»7i </lmu> 2 =»? 2 

Similarly for one has 

£>i(ri = 2u!,xi,r 2 = 2u 2 ,x 2 ) 

1 f f r :E 2-Ar( Mc (A2)-i) 

/ «£i / g£2 ^2 _ _ _ e jy(g^( ft) , /?1 (zi)+ g>te(/3a) , ga (i/ a2 )) 



(2vri) 2 (l - a)2(«i+«2) / ^ J Z2 xi-jv^O-i) 

Zi [\z\ — \ 11 — z 2 

Zl - z 2 1 2*i + 1 2 1 + 2 2 



X _?^Ji*^-II^U, (4.78) 
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where /3» = Ui/N. The critical points of Zi,z 2 are Z\ = p c (/3x),z 2 = l/p c (/3 2 ). To leading 
order, the paths are 



zi = Pc(Po) 1 + 



1 



-=Pc(/3o) 1 
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rf(/3o)iW3 

1 



Tl - ZWi 



r 2 - m> 2/ 



(4.79) 
(4.80) 



so that we have 

D x ~ (1 - a)- 2( " 1+M2) 



-( Pc (/3o))- {6+6Wo)Arl/3 - 



l + p c (A)))(i - p c (A))) 3 v v " ^(/5o) 3 iV 

2ArAo(/3o)+A 1 (/?o)c(/3o)V 2 / 3 (Ti+T2)+A 2 (A>)cW 



4vr 2 



Im«)i=r;i 



dw 2 (n - r 2 - - iu 2 )) e^^+^^+t^^). 



(4.81) 



Imu> 2 =^2 



Now suppose that we substitute the asymptotic expressions (j4.65j) . (j4~73|) . (j4.77]) . (|4.81|) 



into (jH.fijh Notice that the kernel of the form 



a b 
c d 



and 



aa fib 
c//3 d/a 



give the same value 



for the determinant and hence that some of the prefactors in (|4.65 ) 
have no effect on the value of the determinant. Moreover, in our scaling limit, the off- 
diagonal elements, I\ and D\, vanish due to the difference of order in N. The diagonal 
elements are transpose to each other and hence one finds 



F[h(r u 2N) < X x , ■ ■ ■ , h{r m , 2N) < X m ] = y/det(l + K l9 ) 
-> v / (det(l + K, 2 G)) 2 = det(l + K 2 G), 



(4.82) 



where g is g{r h x) = —xj^x) (i = 1,2, ■•• ,m), where Jj = (X^oo) and = 2/3 N + 
2c(p )N 2 / 3 r l , X t = a(/3 + + M^o)^ 3 for i = 1, 2, • • • , m. □ 



4.3 Near the Origin 

In the last subsection, we saw that the fluctuation in the bulk is described by the Airy 
process. On the other hand, as already mentioned, we know that the height fluctuation 
at the origin is given by the GOE Tracy- Widom distribution. In this subsection, we are 
interested in the crossover between them. Let us define the scaled height variable near the 
origin as 

h(r = 2cN 2 / 3 r,t = M = 2N) - aN 2 .„ . 

H N (r) = -i ^ 1 + r 2 , (4.83) 
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where 



a = a(0) 
d = d{0) 
c = c(0) 



2a 
1 — a 

a l/3( 1+Q )l/3 

1 — a 

(1 + a)s 
i 

as 



(4.84) 
(4.85) 
(4.86) 



The second term in f!4.83|) comes from the expansion of a{0) in ()4.27|) around (3 = 0. 

We show that the fluctuation of the model near the origin is described by the process 
which gives the orthogonal-unitary transition in random matrix theory. Namely, we show 



Theorem 4.2. 



lim F[H N {n) < si, • • • , H N (r m ) < s m ] = v / det(l + K£), 



(4.87) 



where G(Tj,£) = — X(s ,oo)(C) (j ' = 1; 2, • • • , m) and K,\ is the 2x2 matrix kernel, for which 
a representative element is given by 



£i(ti,£i;t 2 ,£ 2 ) 



<Si(ti,6i;t" 2 ,6) l»i(ri,Ci;r2,6; 



(4.* 



i/ie matrix elements being 

<Si(t 1; £i; r 2 , £ 2 ) = Sifa, £ l5 r 2 , 6) - $^(6, 6), (4-89) 

5i(ri,ei;r 2 ,6)= / ^Ae- A (— )Ai(^i + A)Ai(6 + A) + -Ai(6) / rfAe^Aife - A), 
Jo 1 Jo 



^i(ti ) 6;t2,6) 



rfAe- Ari Ai(^ - A) / due- UT2 Ai(£ 2 - «) 

J A 



^l(Tl,Ci;7"2,C2) = 



poo /»oo 

/ dAe~ Ar2 Ai(£ 2 — A) / dve^M^i 

JO J A 

1 f 00 d 
-- J dAe-^Ai(6 + A)- {e^Ai^ + A)} 

1 f 00 d 

+ -J dAe-**Ai(£i + A)— {e- A ^Ai(e 2 + A)} . 
Definition of $ T i,t 2 (£i, 62) already given in \4-%4 ^ 



(4.90) 



(4.91) 



(4.92) 



Remark. The same kernel appeared in the context of the orthogonal-unitary transition 
in [2H1 and the problem of vicious walks in |3Sj- For the special case of a fluctuation at 
t = 0, we have 

lim F[H N (0) < s] = Fi(s), (4.93) 

iV— >oo 
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where F±(s) is the GOE Tracy- Widom distribution [22] • It is remarked that this result 
was anticipated in jHj based on the Meixner orthogonal ensemble representation and was 
shown in |25j . On the other hand, for r — > oo, X\ and T>\ goes to zero. We recover the 
Airy process, which is consistent with the results of the previous subsection. 

Proof. Let us derive (POTT)) . We first notice that flOU]) , (|PT|) . P~32Jl have double contour 
integral representation. Assume n, T2 > and ??2 > 0. For «Si, one has 

<Si(ti,£i;t 2 ,£ 2 ) 

(iwi / G?U> 2 I T7 r + 



4vr 2 Ji mwi=ril Ji mW2=m V r 2 - Ti + z(tui + w 2 ) 2(r 2 + z'w 2 ) 
e i£iwi+i6™2+§ (wf+w%)^ (4.94) 

where r/i + 772 + Ti — r 2 > and r 2 — r/ 2 > 0. For Xi, 

^i(n,6;i"2,6) 

' ; d^l / dw 2 - —— — J^l+i^ + ^wf+w^) 



4vr 2 Ji mm=ril Ji mW2 = m fa + iw 2 ) (n + T 2 + i(w!+ w 2 )) 
-(ri,6^r 2 ,6), (4-95) 
where Ti — 771 > 0, r 2 — r\ 2 > 0. For X>i, 

^i(n,6;^,6) 

= JL / dWl / dw 2 - T _±^!i e ^+^ 2 +iW+-D 

47T J Jltn Wl = Vl Jlmw 2 = m 4(Ti + T 2 - l(tWi + W 2 J) 

-(n,6^r 2 ,6), (4.96) 

where there is no additional condition on r 1; r 2 , 771, 771. These are easily obtained by using 
the integral representation of the Airy function, 

Ai(0 = ^ / dwe* w+i ^. (4.97) 

Im w=r] 



2n 

In the following, we show, as iV — > 00, 

5i(r 1 ,x 1 ;r 2 ,x 2 ) ~ e^-^/ 3 ~^ + ^dN 1 / 3 S 1 (r 1 ^ 1 ;T 2 ^ 2 ), (4.98 

h(r 1 ,x 1 ;r 2 ,x 2 ) ~ e-^+^+^+^J^n, ^i; r 2 , 6), (4-99) 

D x {r u x x -r 2 ,x 2 ) ~ e^+^^^-^^iV 2 / 3 ^^, ^; r 2 , 6), (4.100) 

where we set = 2ciy 2/3 Ti, 2; = aiV + (& - r 2 )<iiV 1/3 for z = 1,2. From these, (Q7I) is 
almost obvious. 

Let us first consider the asymptotics of S\. Basically one can follow the same route as 
in the last subsection with /? = 0, but there appear some simplifications and differences. 
For instance, ()4.61|) . (j4.62|) reduce to 

Zl = 1 + divW Tl " lWl) ' v 2 = 1 ~ w^ {T2 + lW2) - (4 - 101) 
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Then we have 



Z\ L, Z\ — Z 2 



dN 1 / 3 1 dN 1 / 3 



1 H > ~ 1 . (4.102) 

zi-z 2 { (1 + z 1 )(z 2 - 1) J r 2 - Ti + i(wi + w 2 ) 2r 2 + iw 2 

Compare this with (|4.64|) . where the second term on the left hand side was negligible. 
Expansion of G(/3) in (J4.58)) is now simply 

(y 

G{(3) = 2 log(l -a) 13+ , 2 /3 3 , (4.103) 

6(1 + a) z 

so that one has 

e M?Mc(0i),/3i (pc (ft)) ~ (l - a) 2Ml e T i 3/3 , (4.104) 

e JVPMc{ft,)./Ja(p»Ofe)) ~ (1 _ a )~ 2u2 e- T 2 /3 , (4.105) 

instead of (|4.59J) . (|4.60J) . Asymptotics of 0, which is already given in [T^], can be obtained 
by setting = in (J4~?5j) : 

n , r2 (x 1? x 2 ) = e^-^^^-^d^ 1 / 3 ^,,,^!,^). (4-106) 

Combining these, we see (|4.98|) holds. 
As for Ii, (|4TF5J) . (|4~77)J) reduce to 

z 1 = l- 1 ^(r 1 + iWl ), 1 = 1-^^ + ^). (4.107) 



Since 



ZL 



lz 2 + l ll + ^il d 2 iV 2 / 3 / -1 -1 



(4.108) 



Zi - z 2 \ 2 z 2 - 1 2 1 — 2! J n + r 2 + z(u>i + iu 2 ) V r 2 + ^2 n + ^i 

we see that (EHIHj) holds. 

Finally, for D u (gZZSD , dUHDD reduce to 

2l = 1+ d7vW (Tl-m) ' i = 1+ div^ (r2 "^ 2) - (4 ' 109) 

Since 

Si f 1 1 ~ z 2 1 _ 1 (n - wi) - (r 2 - w 2 ) f4H0) 

zi - z 2 \ 2 zi + 1 2 1 + z 2 J ~ 4 n + t 2 - i(wi + w 2 ) ' 

we see that (j4.100jl holds. □ 
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5 Model without Source at the Origin 



In this section we study the case where 7 = in (|2.4j) . As already mentioned in the 
introduction, this case corresponds to the PNG model without the source at the origin. 
The analysis proceeds in a fairy parallel way to the orthogonal case. The main difference 
is that there is a new restriction on the height lines hf, h 2i (0,t) = h 2 i + i(0,t) + 1 (i = 
0, 1, • • • , N/2 — 1) at time t = M = 2N. Here and hereafter we assume N is even. This 
condition might look somewhat strange at first sight, but in fact can be understood after 
a short reflection. An example of a multi-layer heights is given in Fig. 3, in which one 
sees that h 2i and /i 2 i+i form a pair at r = 0. As in the orthogonal case, we change 
the way of looking at multi-layer heights. The space coordinate r will be interpreted as 
the time coordinate. The height coordinate will be interpreted as the space coordinate 
and is represented as x. Then the restrictions read xf 1 = 1 — % for % = 1, 2, • • • ,N and 

2& = x|U-lfori = l,2,...,iV/2. 



5.1 Determinantal Process and Kernel 

We consider the weight of non-intersecting paths of n particles given by 

w n , M {x) = det ( M^"!' f j L ) II det (^,r + i«, ^ +1 ))m=i> 

\ 'rO,l\ 2i— 1 ' jJ/ i=l t ... ,n r =l 

where n is assumed to be even and xf 1 (i = 1,2, • • • ,n) is fixed. As in the orthogonal 
case, this is slightly different from the weight of the multi-layer PNG model even if we take 
n = N, but we employ this for the moment and take the infinite particles limit when it 
becomes easy. The partition function and the probability of the non-intersecting paths are 
given by ()3.2|) . (j3.3j) respectively. Again, our main focus is on the special case in f)2.6|) . f)2.7j) . 
but the results of this section do not depend on a specific choice of 0. For instance, when 
particles can hop only to nearest neighbor sites for each time, the determinantal process 
is a variant of the vicious walks with a condition that walkers start in pairs and end at 
neighboring sites |42j . 

Under this arrangement, we can show 

Proposition 5.1. 

M-l n 

E II Hil + g^x^PnMx) = Vdet(l + K,g), (5.2) 

x r=0 j=l 

where g is some function. The determinant on the right hand side is the Fredholm deter- 
minant, where K4 is a 2 x 2 matrix kernel, 



K 4 (r 1 ,x 1 ;r 2 ,x 2 ) 



S i (r 1 ,x 1 ;r 2 ,x 2 ) D 4 (n,x 1 ;r2,x 2 ) 
h(r 1 ,x 1 ;r 2 ,x 2 ) S^r^x^r^xi) 



(5.3) 
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with the matrix elements being 

S 4 (r 1 ,x 1 ;r 2 ,x 2 ) = S 4 (r 1 ,x 1 ;r 2 ,x 2 ) - (f> ri>r2 (x 1 ,x 2 ), (5.4) 

n 

S i (r 1 ,x 1 ;r 2 ,x 2 ) ^ ^ M {x u x^){A^ 1 ) i jG i (r 2 ,x 2 ;M,xf), (5.5) 
J 4 (ri, xi; r 2 , x 2 ) = h{ri,Xi] r 2 , x 2 ) - G^n, x x \ r 2 , x 2 ), (5.6) 

n 

I A {ri,x 1 ;r 2 ,x 2 ) = - ^ G 4 (r x , x l5 M, xf )(Aj 1 )„-G 4 (r 2 , x 2 ; M, xf ), (5.7) 

n 

D 4 (ri,xi;r 2 ,x 2 ) = 0n,A/(xi, xf XAj^y^jv/O^, xf ). (5.8) 



i,3=l 



Here 



(A^) i:j = ^ s (yz ~ 2/i)0o,m(£/i, )0 o ,Af(z/2, xf ), (5.9) 

2/1,2/2 

G 4 (r 1 ,xi;r 2 ,x 2 ) = s(y 2 - yi)4>o,n{Vi: ^i)<?V 2 (y2, x 2 ), (5.10) 



2/1,2/2 



'J 



s(« - j) = - ( 5 - 1] 



Remark. The subscript 4 refers to the fact that this case is related to the symplectic- 
unitary transition in random matrix theory. It should also be remarked that at this stage 
the infinite particles limit is taken easily. One only replaces the summation, Y^ij=n m 

(E3.0).f|lHi it ii J2Z=v 

Proof. As in the orthogonal case, we derive ()5.2|) by generalizing the methods of [T9*l l34j . 
Using the Heine identity (|3.16j) . the partition function Z n ^\g\ is written as 

^bi^E^K^riS-SM ■ (5 - i2 » 



«=!,■■■ ,n/2, j=l,- - ,n 



Here 



<«/(•' ''--''j') = J2 0- + 9(0^))Mx° i ,X 1 )(l + g(l,X 1 )) ■ ■ -0M-i,M(*M-i,xj 



My 
b 3 

Xi,— ,Xm—i 

= o ,A/(x°,xf ) + (0 o ,n • g^){xlxf). (5.13) 

In the second equality we have used ipj defined in (I3.28J1 . From (|5.12j) . ()5.13|) and the 
identity 

v 

:{<l>i(Xj)i)i(Xj))i=i...2n,j=l---n I 

\ 2n 

((2n)!) 2 det ( ^(Hv^jiv) ~ <t>jWi(v)) J , (5-14) 

*,i=i 



( £ det < 

\a;i,— ,x n 
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one finds 



where A 4 is defined in ()5.9|) and 

(4 1} ) M - = ^[(0 o ,n ■9il> j )(y,x?)<t> Q , M (y + l^f) -MVtXkKfarr -9^)(y+ hxf)} 
y 

= Y J gA-(G A ) 3 (5.16) 

(4 2) ) M - = ^[0 o , A f(i/,xf )(0o,n ■ + ) - (<Kn ■ gip k )(y,xF)<l>oMv + )1 

2/ 

= -X)^--(G 4 )i (5.17) 



|(3)> 



(5.15) 



- (0o,n -g^kjiy^k )(0o,n •0Vj)(y + 1 >zJ tf )] 

- ' G 4(#i)- 



(5.18) 



We should notice that the forms of , A 4 2 ^ and A^ 3 "* are the same as the corresponding 
ones in the orthogonal case, ()3.29j) . (|3.31j) and ()3.32j) . respectively. Only difference is that 
sgn in (jHH2| and ()3.13j) is replaced by s, (j5.11|l . in (|5.9jl and (|5. !())) ■ Thus we can calculate 
the kernel along the line in section 3. The result is 



Z n ,M[g] 

= det | 1 + 



Recalling the definitions of 0j and (G^)i, we see that (|5.2jl holds. 

5.2 Scaling Limit 

5.2.1 Generating Functions 

The symbol a 4 (z) of the matrix v4 4 is computed as 



(5.19) 
□ 



a*{ z ) = fo,M - s 4 (z)f 0;M {z) 



where 



s 4 (z) = ^2 s(k)z k = -z + -. 



(5.20) 
(5.21) 



feez 
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The winding number of a±{z) is again not zero, but the difficulty can be overcome as in 
the orthogonal case, s^z) can be expressed in two ways 

s,{z) = s ( ll(z)s^(z) = sfl(z)sfl(z), (5.22) 

with 

s$(z) = l + z, s£(z) = -l + -, (5.23) 

z 



sfl(z) = l-z, s fl(z) = l + -. (5.24) 

z 

Then in terms of a^, a!£L for i — 1, 2 defined by 

a£(z) = /o,m ; - (£) (5-25) 

at(z) = / 0iM;+ / ,m;-(^2(*), (5-26) 

1 is expressed as 

^^{ T fe) r fe) +T (i) T (i ir (527) 

Using the generating functions of 1 and G 4 , 



,3 Z 2 



Z 1 1 1 



* - * w(^i) w(^) w(^)/o^ ; -(i) 2 \ 4^)4^) 

(G , 4 ) n , r2 (^) = -/o,ri(^)s4(^)/o,r 2 Q , (5-28) 
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we can calculate the generating functions of £4, J4, D4; 

§4(n, 21; 7*2, 22) = ^(^^i!^^)^ 1 ;^ 2 



«1 



fr U MAi)krrABkM-A^) i f 45(ir) , 42(i) 



*1 - ^2 /0,M;-(^l)/0,ri; + (^)/r 2l M i -(i) 2 \ s «(l) + s g(| 
(1 _ q)2(u 3 -u 1 )( 1 _ q/ Zl )^+»l(l - aZ /^ ^ f ^-Z 2 { 

(l-a Zl )^(l-a/z 2 )N+u 2 z 1 -z 2 \ + z 2 {l-zl)<' ' 

I 4 (r 1 ,z 1 ;r 2 ,z 2 ) = ^ h(n, Xi,r 2 , x 2 )z^ z^ 2 

Z\ /o,r i; +(^l)/o,M;-(j-)/o,r 2 ;+(^)/o,A/;-(^2) 



(l-a) 2 ^+^{l-a/ Zl ) N -^{l-az 2 ) N -^ z x f 1, , r 

z\-\ > , (o.ouj 



{l-az 1 ) N + u i{l-a/z 2 ) N+u 2 z 1 -z 2 { z 2 

^4(7*1, zr, r 2 , -83) = ^4(7*1,^1; 7*2, ar 2 )^ 



2l 



£1 ^—£2 

fr 1) M;-('t)fra,M;-( z ^) if 1 1 



+ 



*i - -2 i^A^S^Ak^A^U^A^ 2 I sfli^sfK^ 
{l-a/z 1 ) N+u i{l-az 2 ) N+u * 



(1 - a) 2 («i+^)(i _ azx) N - u ^{l - a/za)^"" 2 



— z 2 y 2 2Tx + 1 1 — Z 2 21 — Zxl + z 2 
5.2.2 Bulk 

In the bulk region, the asymptotic behavior is just the same as in the orthogonal case. 

Proposition 5.2. For the scaled height variable defined by \4.2b\ l, M-SUj) holds in the 
symplectic case as well. 

Proof. To show this, we need the asymptotics of S 4 , 7 4 and -D 4 in the bulk region. Applying 
the same reasoning as in 4.2, one finds 

S A ~ (1 - «) 2( -- w ^P c (^o)) K2 ^ l) ^^ 1/3 ^^ 

e Gi(/3 )c(/3 )iVV3( Tl _ r2 ) + ^ 

At / d Wl I dw 2 ( \ ^ e »a»i+»^2+i(»?+^). ( 5 . 32 ) 

4tt 2 Jj mwi=r]1 Jimw 2 =m V r 2 - n + + w 2 ) J 
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4 1 ] (l-Pc(-f3 ) 2 y Pc{P0>> d((3 yN 

e -2WA (A>)-Al(ft)c(ft)i^ /3 (n4^ 

\l dwi I dw 2 (n - r 2 + - iu 2 )) e ^+^+i«+™I) , (5.33) 

»/Im«ji— m Jlmw}.=m 



4tt 2 



e 2ArA (AO+MA))c(/3o)A r2/3 (Ti+T2)+A 2 (A))^ 
1 

4rr 2 

Then discussions similar to ones below (j4.81|) lead to (j4.3()jl . □ 



/ dwt / dw 2 (n - r 2 - - tw 2 )) e ^+^+f K+™f). (5 34) 

' Imuii =rji •/ Im«)2=12 



5.2.3 Near the origin 

Next we show that the fluctuation near the origin is described by the process which gives 
the symplectic- unitary transition in random matrix theory [21] • Namely, if we define the 
scaled height variable as (|4.83|) . one can show 

Theorem 5.3. 



lim F[H N {n) < s u ■ ■ ■ , H N {r m ) < s rn ] = v / det(l + /C 4 £), (5.35) 

N— >oo 

where G{Tj,£) = — X(sj,oo){Q (j ' — 1 3 2, • • ■ , m) and /C 4 is the 2x2 matrix kernel, for which 
a representative element is given by 

<S 4 (ti,£i;t 2 ,£ 2 ) ^(n, £1; T2, £2) 

Z4(ti,6l;t2,6) ^4(7-2, 6; n, 6) 



(5.36) 



^4(Tl,6; ^2,6) 

i/ie matrix elements being 
5 4 (ri,^i;r 2 ,6) = S 4 (r 1; 6; 7" 2 , 6) - $n,T 2 (£i, 6), (5-37) 

/•OO 1 /*oo 

4(^,6^2,6)= / d\e~ x ^ )Ai(6 + A)Ai(6 + A) - -Ai(£ 2 ) / dAe _ATl Ai(£ 2 - A), 

./o ^ J£i 

(5.38) 

/•OO J 

X 4 (r 1 ,6;r 2 ,6) = - J dXe-^Mfo - A)— {e^Aife - A)} 

r°° d 

+ J dAe^Aife-A)— (e^Ai^-A)} (5.39) 

1 /*oo /*oo 

^4(ri,6;r 2 ,6) = 5 y o rfAe- AT1 Ai(6 + A)^ ^e^Aife + u) 

i /*oo 

dAe- Ar2 Ai(^ 2 + A) / dvAi(£i+u). (5.40) 



1 

4 Jo 



A 



Definition of $ n ,r 2 (Ci) 62) already given in 
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Remark. For r = 0, we have 



lim F[H N (0) <s} = F 4 (s), (5.41) 

N— >oo 



where F^s) is the GSE Tracy- Widom distribution [22] • Notice a notational difference 
in [22] and [2S]- We follow the convention in the latter; our F±(s) is F^s) = Ff R (s) = 
F™(V2s). The result, EQ) , was shown in [2E]; it would also be possible to prove this 
by using 



JV 

_> 



F[h(0,2N)<X] = -@ Yl II {h i -h j ) 2 (h i -h j + l)(h i -h j -l)]lq h *, 

Zn h6N * !<^<f 1=1 

max{hj}<X+N-l 

(5.42) 

and the skew orthogonal polynomials techniques [ISIEIE^]- The Meixner symplectic en- 
semble representation, ([5.42)1 . was not given in [B] but can be proved similarly if one notices 
that a symmetric N x N matrix with zero elements on diagonal and non-negative integer 
elements on off-diagonal is mapped to a semistandard Young tableau with all columns of 
even length through Knuth correspondence [IB]. 

Proof. Derivation of ()5.35[) is analogous to that of (|4.87j) . Using the integral representation 
of the Airy function, ()4.97j) . we have, for «S 4 , 

5 4 (ri,^i;r 2 ,5) 

^ dwi [ dw2 



4?r 2 Ji mWl=m Jimw 2 = m V r 2 - n + i(w! + w 2 ) 2(n - iwi) 

where + ?7 2 + 7i — r 2 > 0. For 2T 4 , 
Xi(ti,Ci;t2,6) 

1 /" dwi f d Tl + iw2 to +^ 2 +f(-?+»§) 



4tt 2 Ji mwi = m Jiu, W2=m On + r 2 + i(wi + w 2 )) 
-(ti,6^t 2 ,6), (5-44) 

where n + r 2 — 771 — r/ 2 > 0. For V4, 



2 



' 1 dw X [ dlU 2 - — J^+i^ + ^wf+wl) 



4tt 2 J lmwi=m Jim W2 = m 4(n - iwi)(ti + r 2 - f(^i + w 2 )) 
-(ti,£i<->t 2 ,6), (5-45) 



34 



where there is no additional condition on ti,T2,T]i,T]i. Then applying the same method as 
in 4.3, we can show 

S 4 (r 1 ,x 1 ;r 2 ,x 2 ) ~ e^-^+^-^dN^S^^r,,^), (5.46) 
h{r 1: x lV r 2 ,x 2 ) ~ eW^J/K^v^n.^T^), (5.47) 
L> 4 (ri,x i; r 2 ,x 2 ) ~ e-^ 3 )/ 3+ ^ ri +^P 4 (Ti,6;T2,6), (5.48) 

from which ()5.35|) follows. □ 

6 Continuous Limit 

Let us take the limit a — > 0, A" — > oo with i = aA^, fixed. In this limit time and space 
become continuous; the time variable t can take any positive value whereas the scaled space 
variable is defined to be v = ar/2. The model is reduced to the standard PNG model in a 
half space (v > 0) with an external source at the origin. As in the discrete case, we extend 
the space to the whole space (u 6 1), putting the symmetry condition on the hight with 
respect to the origin. First we have a flat substrate. At time 0, a nucleation of a height one 
occurs at the origin. This step grows laterally in both directions with unit speed. Above 
this ground layer there occur other nucleations with rate two per unit length. The height 
of a nucleation is always one. There is an external source at the origin with rate 7. As in 
the discrete case, one can define the multi-layer version [T3j . 

It is easy to take this limit at the level of the generating functions. For the orthogonal 
case, (|I22l,{£23l,{£23|l become 

§ x rv, c 2(vi-V2) c t(zi-l/z 1 )-v 1 (z 1 +l/z 1 ) c t(l/z2-z 2 )+V2(z2+l/z 2 ) Z l L + Z l ~ Z 2 

z\ - z 2 \ (1 + z 1 )(z 2 - 1) 



(6.1) 

| rs. c -2K+^2)^i-lM)+fi(^i+lAi) r .f(l/^-^)+^2^2+lA2) Zl \lZ 2 + l 11 + Zl 



z\ — z 2 \2z 2 — l 21 — zi 



D ^ e 2(vi+V 2 ) e t{zi-l/zi)+Vl(zi+l/zi) e t(l/Z2-Z2)+V2(z2+l/z 2 ) Zl I' 1 " ! ' ' '~- 



zi - z 2 I 2 zi + 1 21 + z 2 



(6.2) 



(6.3) 



Expressions for the symplectic case are given by (j6.1jl . (|6.2jl . (|6.3jl with the last factors 
replaced by the corresponding ones in (j5.29j) . (j5.3U|) . (j5.31|) . 

The scaling limit can also be studied as in the discrete PNG model. First the thermo- 
dynamic shape is known to be 



h(v = o t,t)/t~2yJl-0*, (6.4) 
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where < /3q < 1 is fixed pQ. To consider the scaling limit in the bulk, define the scaled 
height variable as 

H n (t,(3 ) = — — -i l - , (6.5) 

d(J3o)t3 

where 



a{(3) = 2y/l -/3 2 , (6.6) 

d(P) = (1 - p)l , (6.7) 

c((3) = (1 - P 2 )l (6.8) 

Then one can show 

Proposition 6.1. For the scaled height variable defined by 16'. 5\) . M-HUp holds. 

The proof of the Proposition 16.11 is parallel to that in 14.21 The main difference is that 
the function g^p in (|4.42j) is replaced by a simpler 

1 1 

9nA z ) =z + -) - filogz, (6.9) 

z z 

for which the double critical point is 

(i c (j3) = 2y/T=]P, (6.10) 

PM = J\^-- (6.11) 

As for the fluctuation near the origin, if we define the scaled height variable as 

. , h(v = £sr, t) — It 9 . 
H n (t) = ^ + r 2 , (6.12) 

U 

we have 

Proposition 6.2. For the scaled height variable defined by one has li4-87]) for the 

orthogonal case and \5. 35]) for the symplectic case. 



7 Discussions 

In the preceding sections, we have given a detailed analysis of the height fluctuation of 
the model for two special values of the strength of the external source at the origin, 7 = 1 
and 7 = 0. Our results show that the fluctuation near the origin is described by the 
orthogonal/symplectic to unitary transition in random matrix theory. 
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This implies, in particular, that the height fluctuation of the PNG model at a single 
point near the origin is equivalent to that of the largest eigenvalue of the transition en- 
semble. To check this, we performed Monte-Carlo simulations of the PNG model and the 
transition random matrix. In Figs. 4 and 5, we have shown the fluctuations of the height of 
the PNG model and the largest eigenvalue of the transition ensemble. We see an excellent 
agreement between them. 

Unfortunately properties of the transition ensemble have not been well studied com- 
pared to the ensembles with a specified symmetry. For example, it seems difficult to 
numerically compute the probability distribution function of the largest eigenvalue . This 
is sharply contrasted to the situation for the GUE/GOE/GSE, for which the Painleve rep- 
resentation allows us to plot the probability distribution function with high accuracy. It 
would be desirable to study the transition ensemble to better understand the statistical 
properties of the PNG model. 

Now we argue what happens for other values of 7. For the fluctuation at the origin, 
the results have been obtained in where a strong universality is observed. The GSE 
Tracy- Widom distribution describes the fluctuation not only for 7 = but also for all 
values in < 7 < 1. The fluctuation becomes the Gaussian for all values in 7 > 1. On 
the other hand, the GOE Tracy- Widom distribution appears only at 7 = 1. The Gaussian 
fluctuation for 7 > 1 is stated only for the continuous model in |2H], but is expected to 
persist for the discrete model as well. In Fig. 6, typical shapes of the droplet are shown 
for three cases where 7<1,7 = 1,7>1. From these one should be able to guess that 
the shape looks similar for all values in 7 < 1 but a cusp-like piece appears near the origin 
when 7 becomes greater than unity. 

This suggests that the fluctuation near the origin also has a similar universality. For 
all values in < 7 < 1, it is expected to be described by the symplectic-unitary transition. 
This is a conjecture because we do not know how to prove this at present, but is supported 
by a Monte-Carlo simulation. See Fig. 7, where a good agreement is observed. 

For 7 > 1, we can also expect some universal behavior for the fluctuation near the 
origin. We have not found a compact formula for general 7, but the situation becomes 
quite simple in the limiting case where a — > with < 7a < 1 fixed. In this limit, the 
nucleations in the bulk are so rare that the only those at the origin are important. Then 
the height at a position r and at time t would be almost the same as that at a position 
and at time t — r. Therefore the height fluctuation at a single point is given by the 
Gaussian. If we set r = pt with < p < 1 fixed, we have 

P ( Mr-O-d-rtM < \ = 1 f e _ 4d5 = 

where ac = jz^ and do = ^ ■ A result of a Monte-Carlo simulation corresponding to 
this case is shown in Fig. 8. In addition, since nucleations at the origin are independent 
for each time, the multi-point equal time height fluctuation would be described by the 
one-dimensional Brownian motion. The situation would be somewhat more difficult for a 
smaller 7, but we expect that the same Gaussian fluctuation is observed in an appropriate 
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limit. 



8 Conclusion 

In this paper, we have studied fluctuation properties of the one-dimensional polynuclear 
growth (PNG) model in a half space with an external source at the origin. We have mainly 
considered the model in a discrete space and time, but have also given results for the model 
in a continuous setting. The results in the scaling limit are the same for both cases. 

For two special values of the strength of the external source, 7, we have performed 
a detailed analysis. The 7 = 1 case corresponds to a critical point, which we call the 
orthogonal case. The 7 = case corresponds to the model without the external source, 
which we call the symplectic case. The main results are f)4 . 8 7|) for the orthogonal case, and 
(|5.35J) for the symplectic case. According to these the height fluctuation of the model near 
the origin is equivalent to those of the largest eigenvalue of the orthogonal/symplectic to 
unitary transition ensemble at soft edge in random matrix theory. We have also shown 
that the height fluctuation at bulk is described by the Airy process. For other values of 
7, we have conjectured that the fluctuation is the symplectic- unitary type for < 7 < 1, 
whereas it is the Gaussian type for 7 > 1. Some Monte-Carlo simulation results are also 
presented to confirm our results and conjectures. 
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Figure Captions 



Fig. 1: Dynamical rules of the discrete PNG model, (a) At time t, a nucleation of a height 
k occurs at site r with probability (1 — q)q k when r ^ and (1 — 7 A /^)(7v / ^) fc when r = 0. 
Once the nucleation occurs, it grows laterally toward right and left with unit speed, (b) 
When two steps collide, the one with higher height swallow the one with lower height. In 
the multi-layer version, a nucleation occurs in the lower layer at the position above which 
the collision occurs in the upper layer. The height is equal to the height of the swallowed 
region in the upper layer. 

Fig. 2: An example of the multi-layer PNG heights for the 7 = 1 case for q = 0.25 at (a) 
t = 6 (b) t = 100. It looks symmetric with respect to r = 1/2. The limiting shape is also 
shown clS db dotted line for (b). 

Fig. 3: An example of the multi-layer PNG heights for the 7 = case for q = 0.25 at 
t = 100. Each neighboring two heights form a pair at r = 0. The limiting shape is also 
shown dotted line. 

Fig. 4: Probability distributions of the scaled height of the PNG model for 7 = 1 and 
the scaled largest eigenvalue of the orthogonal- unitary transition random matrix ensemble. 
For the PNG model, the parameters are for q = 0.25, t = 2000, and (a) r = (b) r = 330 
(c) r = 1000 (each 10000 samples). The first two cases correspond to (a) r = (b) r = 1, 
in which the scaling ()4.83|) is used. The case (c) is better fitted to the bulk scaling limit; 
the scaling ()4.26|) is used. In the figures, they are represented as circles. For the transition 
random matrix ensemble, the data are for N = 500, (a) r = (b) r = 1 (c) r = 10 (each 
1000 samples). In the figures, they are represented as +. For (a)(resp. (c)), the GOE 
(resp. GUE) Tracy- Widom distribution is also shown as a solid line. 

Fig. 5: Probability distributions of the scaled height of the PNG model for 7 = and 
the scaled largest eigenvalue of the symplectic-unitary transition random matrix ensemble. 
For the PNG model, the parameters are for q = 0.25, t = 2000, and (a) r = (b) r = 330 
(c) r = 1000 (each 10000 samples). The first two cases correspond to (a) r = (b) r = 1, 
in which the scaling ()4.83|) is used. The case (c) is better fitted to the bulk scaling limit; 
the scaling (|4.26J) is used. In the figures, they are represented as circles. For the transition 
random matrix ensemble, the data are for N = 200, (a) r = (b) r = 1 (c) r = 10 (each 
10000 samples). In the figures, they are represented as +. For (a) (resp. (c)), the GSE 
(resp. GUE) Tracy- Widom distribution is also shown as a solid line. 
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Fig. 6: Typical droplet shapes for three values of 7 = 0.8, 1.0 and 1.2, which are fairly 
close to the critical value 7 = 1. The model parameters are taken to be t = 1000, q = 0.25. 
Notice that the appearance of the shape changes drastically when 7 gets greater than one. 

Fig. 7: Probability distributions of the scaled height of the PNG model for 7 = 0.5 and 
the scaled largest eigenvalue of the symplectic-unitary transition random matrix ensemble. 
Parameters are the same as those for Fig. 5. 

Fig. 8: Probability distributions of the scaled height of the PNG model for q = 0.0001, 7 = 
50 and the Gaussian distribution (solid line). For the PNG model, the parameters are 
t = 2000 and r = (+), r = 330 (o), r = 1000 (•) (each 10000 samples). 
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